Supercurrent- Induced Magnetization Dynamics 
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We investigate supercurrent-induced magnetization dynamics in a Josephson junction with two misaligned 
ferromagnetic layers, and demonstrate a variety of effects by solving numerically the Landau-Lifshitz-Gilbert 
equation. In particular, we demonstrate the possibility to obtain supercurrent-induced magnetization switch- 
ing for an experimentally feasible set of parameters, and clarify the favorable condition for the realization of 
magnetization reversal. These results constitute a superconducting analogue to conventional current-induced 
magnetization dynamics and indicate how spin-triplet supercurrents may be utilized for practical purposes in 
spintronics. 
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Introduction. The interplay between superconducting and 
ferromagnetic order is presently generating much interest in 
a variety of research communities [1 1. Besides the obvious 
interest from a fundamental physics point of view, a major 
part of the allure of superconductorjferromagnet (S|F) hy- 
brid structures is the prospect of combining the spin-polarized 
charge carriers present in ferromagnets with the dissipation- 
less flow of a current offered by the superconducting envi- 
ronment. By tailoring the desired properties of a hybrid S|F 
system on a nanometer scale, this interplay opens up new per- 
spectives within spin-polarized transport. 

Closely related to the transport of spin is the phenomenon 
of current-induced magnetization dynamics Q. The gen- 
eral principle is that a spin-polarized cuiTent injected into a 
ferromagnetic layer can act upon the magnetization of that 
layer via a torque and thus induce magnetization dynamics 
IJI. Previous works in this field have considered mainly non- 
equilibrium spin accumulation via quasiparticle spin-injection 
from a ferromagnet into a superconductor |l4]-(8l . The concept 
of supercurrent-induced magnetization dynamics suggests an 
interesting venue for combining the seemingly disparate fields 
of superconductivity and spintronics. However, magnetiza- 
tion dynamics in a Josephson junction has so far been dis- 
cussed only in a handful of works ll9 HT3l . In particular, it 
was demonstrated in Ref. Q how a supercurrent would in- 
duce an equilibrium exchange interaction between two non- 
coUinear ferromagnets in an S|F|N|F|S junction (N stands for 
normal metal). Taking into account the fact that a Joseph- 
son current flowing through such an inhomogeneous magne- 
tization profile will have a spin- triplet contribution UJ, such 
an interaction implies that it should be possible to generate 
supercurrent-induced magnetization dynamics in this type of 
junction. This would constitute a superconducting analogue 
to magnetization dynamics in a conventional spin-valve setup. 
To this date, this remains unexplored in the literature. 

Motivated by this, we study in this Letter for the first time 
the magnetization dynamics of a multilayer ferromagnetic 
Josephson junction when a cuiTent bias is applied, based on 
the Landau-Lifshitz-Gilbert (LLG) equation |14|. Our main 




FIG. 1: (Color online) The proposed experimental structure. Two 
ferromagnetic layers separated by a normal spacer are sandwiched 
between two i-wave superconductors. The magnetization directions 
may be non-aligned and tuned by means of an external field, as long 
as the exchange coupling between the ferromagnets is sufiiciently re- 
duced by the thickness of the normal spacer. By current-biasing this 
system one may generate a supercurrent which induces magnetiza- 
tion dynamics. 



idea is to utilize the spin-triplet nature of the Josephson cur- 
rent, which very recently has been observed experimentally 
ifTSl - EOl , in order to induce a torque on the magnetic order pa- 
rameter and thus generate magnetization dynamics. The ex- 
peiimentally relevant setup is shown in Fig. [T] two feiTomag- 
nets with different coersive fields are separated by a normal 
spacer and sandwiched between two conventional ,s-wave su- 
perconductors. The coersive fields are such that the magnetic 
order parameter is hard in one layer, while soft in the other 
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FIG. 2: (Color online) The time-evolution of the normalized magnetization components nij. Here, we have set ^ = 290 with an initial angle 
misalignment Bq/ti = 0. 1. We have considered the weak damping regime a <C 1 in the upper row, setting (a) CO = 0.5 and a = 0.05 while in 
(b) CO = 1 .5 and a = 0.05. In the lower row, we have considered stronger Gilbert damping and set (c) O) = 0.5 and a = 0.5 while in (d) CO = 1.5 
and a = 0.5. The insets display a parametric plot of the y- and z-components of the magnzetization with the red circles indicating the point 
t = (the initial magnetization configuration). The time-span in insets is T e [0, 50]. 



By application of an external field Hext, it is thus possible to 
tune the relative orientation of the local magnetization fields in 
the two layers. When this junction is current-biased, a super- 
current flows without resistance up to a critical strength. The 
supercurrent strongly modifies the equilibrium exchange in- 
teraction between the ferromagnetic layers and should thus be 
expected to result in supercurrent-induced magnetization dy- 
namics. We investigate this by solving numerically the LLG 
equation, which provides the time-evolution of the magneti- 
zation components in the soft magnetic layer. As we shall 
see, a number of interesting opportunities arise in terms of the 
torque exerted on the soft layer by the Josephson current. We 
proceed by first establishing the theoretical framework used in 
the forthcoming analysis and then present our main results. 

Theory. The magnetization dynamics of the Josephson 
junction is governed by the LLG equation, which in the free 



magnetic layer reads: 
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where y is the gyromagnetic ratio and a is the damping con- 
stant. The effective magnetic field is obtained from the free 
energy functional via the relation Heff = ~ (S^/Swil) / ( I^Mq) 
where 1^ is the unit volume and Mq is the magnitude of the 
magnetization. Note that only transverse magnetization dy- 
namics are described by the LLG equation due to the cross- 
terms coupling to TOi, meaning that \mi\ = 1 is constant in 
time. To study the magnetization dynamics induced by the 
supercurrent, we must identify the effective field. To this end, 
we employ the following phenomenological form |9| stem- 
ming from the exchange interaction generated by the Joseph- 



son current: 
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where A is the unit area, Xf is the Fermi wavelength, and Aq 
is the gap magnitude. The parameters J7, are analogues to 
the quadratic and biquadratic coupling constants for a mag- 
netic exchange interaction. In addition to Eq. (J2]), one should 
also include the anisotropy contribution Fm to the free energy 
which provides an effective field H^ — (Kniy / M^jy where 
K is the anisotropy constant. We have here assumed that the 
anisotropy axis || y. After some algebra, one arrives at the final 
form of the LLG equation in our system: 
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where (|)(f ) = coy? and d is the thickness of the ferromagnetic 
layer. Importantly, we note that the anisotropy contribution 
proportional to K was not included in the effective field used 
in previous works |9|. This contribution is nevertheless essen- 
tial since the supercurrent-induced torque must overcome the 
anisotropy contribution in order to switch the magnetization 
orientation. In what follows, we will present a full numeri- 
cal solution of this equation to investigate the supercurrent- 
induced magnetization dynamics. To this end, we first estab- 
lish experimentally relevant values of the parameters in Eq. 
p). The above equation may be cast into a dimensionless 
form by introducing cOf = jK/Mq, 1 = (apt, CO — (i>j/(i>F, 
and L, = Aq/ {KdXjf). Here, cOf is the ferromagnetic reso- 
nance frequency and (Oy is the Josephson frequency. Employ- 
ing a realistic estimate lfT6l for transmission probabilities in 
the F|N|F part of the system, one finds |9| that Eq. (J2| ac- 
counts well for the Josephson current when j?i = 0.007 and 
J/2 = 0.025. For a permalloy with weak anisotropy, one may 
estimate [JSJ /T ~ 4 x 10"^ K A""*. Moreover, we set Aq = 1 
meV, A,F = 1 A, and d = 10 nm as standard values ffJ] for 
the hybrid structure under consideration. The Josephson fre- 
quency coj is typically of order GHz, but may be tuned ex- 
perimentally. We will therefore consider several choices of 
CO = COy/cO/r and the damping constant a to model a variety of 
experimentally accessible scenarios. As long as the require- 
ment h(Oj <ti Tc is fulfilled, i.e. the time-dependent part is a 
small perturbation, one may consider (|)(f) as a time-dependent 
external potential in the static expression for the free energy 
ifTsl . Since typically coy ^ 1 /jV, this condition is easily met. 
The supercurrent-induced magnetization dynamics come into 
play when the local magnetizations {mi,mR} are misaligned 
with an angle 9o = 9(f = 0) to begin with, and vanishes when 
60 = {0,7t}. The magnetization in the right (hard) layer is 
fixed at niR — y. 

Results and Discussion. Using the parameters discussed 
above, we now proceed to investigate the resulting magneti- 
zation dynamics for both weak damping (a ^ 1) and more 



considerable damping (a ^ 1). Choosing a small initial an- 
gle of misalignment Qq/k = 0.1, the numerical solution of 
the LLG-equation is shown in Fig. |2] The upper row shows 
the weak-damping regime for two different choices of the fre- 
quency ratio (CO < 1 and CO > 1). The qualitative behavior is 
similar in (a) and (b): the supercurrent-induced torque induces 
oscillations of the magnetization components, but is unable to 
switch the magnetization direction in the soft layer from its 
original configuration m^ ~ y. The oscillations eventually die 
out and the magnetization orientation of the two layers sat- 
urates in a parallel configuration. This happens on a faster 
time-scale in (b) due to the larger value of CO used compared 
to (a). Consider now the case of stronger Gilbert damping 
shown in the lower row of Fig. |2] The a.c. Josephson current- 
induced torque now induces more "violent" oscillations and 
is able to rotate the magnetization orientation by K. When 
the frequency is large enough, exemplified in (d) by CO = 1.5, 
full magnetization reversal is achieved and maintained. On 
a larger time-scale, the same behavior is seen for lower fre- 
quencies. In general, we observe in our numerical simula- 
tions that the main difference between the weak and strong 
Gilbert damping regime is that the magnetization dynamics is 
periodic and oscillating in the former case, whereas it tends 
to saturate to a fixed orientation in the latter case. This is 
reasonable as the damping effectively leads to a more rapid 
decay of the oscillating magnetization dynamics induced by 
the a.c. current. The results in Fig. l2] suggest that it is pos- 
sible to generate supercurrent-induced magnetization rever- 
sal in the setup shown in Fig. [T] This phenomenon should 
be intimately linked with the spin-triplet correlations present 
in the non-coUinear setup considered here since spin-singlet 
correlations cannot carry torque. Recently, controllable long- 
range triplet supercurrents have been experimentally observed 
in Josephson junctions with strong inhomogeneous ferromag- 
netic layers I.18-. - .20J . In the present setup, torque carried by 
the supercurrent does not become long-ranged ||2TI (an addi- 
tional source of triplet correlations would be required for that 
purpose |22|) whereas it is still due to spin-triplet, thus me- 
diating magnetic correlations between the two ferromagnetic 
layers. 

We proceed by investigating the role of the anisotropy pref- 
erence in the soft magnetic layer. This is modelled by the 
term proportional to K in Eq. ^. The parameter C, effectively 
models the relative weight of the anisotropy energy and the 
Josephson energy related to the presence of superconductiv- 
ity. In Fig. [3? a), we focus on the behavior of the magnetiza- 
tion component niy. We plot its time-evolution for increasing 
values of ^. It is seen that supercurrent-induced magnetiza- 
tion reversal occurs for ^ > ^^ with a critical value ^j. (for 
the present parameters, we find ^^ — 163). The magnetization 
then saturates, and the magnitude of the oscillations decreases 
as C, becomes larger. In Fig. Brb)-(d), we give a phase-diagram 
for the magnetization switching by plotting my{t -> 00) in the 
a — ^, C ^ CO, and a — CO plane. In general, the results indi- 
cate that the torque generated by the a.c. Josephson effect 
can reverse the magnetization orientation when the Gilbert 
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FIG. 3: (Color online) (a) The time-evolution of the normalized mag- 
netization component my. Here, we have set CO = 1.1 and a = 0.4 
with an initial angle misalignment 9o/it = 0.1. Above a critical 
value of ^, permanent switching occurs from mi \\ y to mi \\ {—y). 
The phase-diagram for magnetization reversal [from the parallel (P) 
configuration mi \\ m^ to the anti-parallel (AP) configuration mi \\ 
(— m«)] is shown in (b-d). A contour-plot is given for m,,(/ -^ <») 
with (b) CO = 0.5, (c) a = 0.5, and (d) ^ = 200. 



damping is non-neglible, the anisotropy contribution is suf- 
ficiently weak and Josephson frequency is sufficiently small. 
The phase-diagram in (b) and (d) shows that the magnetization 
can be trapped in the switched direction when the damping a 
is sufficiently large, in agreement with the results in Fig. 2. 
It should also be noted that the switching becomes less viable 
as CO increases as seen from (c) and (d). When cOy ^ COf , the 
torque from the Josephson current oscillates very fast com- 
pared to the characteristic motion of the magnetization. In this 
way, the magnetization would experience an averaged torque 
over many oscillations, which results in small effect due to 
partial cancellation of the net torque. 



Conclusion. In summary, we have investigated 

supercurrent-induced magnetization dynamics in a Josephson 
junction with two misaligned ferromagnetic layers by solving 
numerically the Landau-Lifshitz-Gilbert equation. We have 
demonstrated the possibility to obtain supercuiTent-induced 
magnetization switching for an experimentally realistic 
parameter set. It is clarified that for the realization of mag- 
netization reversal, large Gilbert damping, small anisotropy 
contribution or a small Josephson frequency are in general fa- 
vorable. These results constitute a superconducting analogue 
to conventional current-induced magnetization dynamics and 
indicate how spin-triplet supercurrents may be utilized for 
practical purposes in spintronics. 
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